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Abstract 

By means of a diagram technique for Hubbard operators we show the ex- 
istence of a spin-dependent renormahzation of the locahzed levels in an in- 
teracting region, e.g. quantum dot, modeled by the Anderson Hamiltonian 
with two conduction bands. It is shown that the renormalization of the levels 
with a given spin direction is due to kinematic interactions with the con- 
duction sub-bands of the opposite spin. The consequence of this dressing of 
the localized levels is a drastically decreased tunneling current for ferromag- 
netically ordered leads compared to that of paramagnetically ordered leads. 
Furthermore, the studied system shows a spin-dependent resonant tunneling 
behaviour for ferromagnetically ordered leads. 

PACS numbers: 72.25.-b, 73.23.-b, 73.63.-b 

Spin-dependent tunneling |I[ and tunneling magneto resistance (TMR) [g,^ have re- 
cently been studied extensively. Concerning spin-dependent tunneling through a quantum 
dot (QD), or similar interacting regions, the main focus has been to investigate the effects of 
a magnetic field applied over the interacting region The opportunity of changing the 

magnetic properties of the leads, leading in and out of the QD, by an external magnetic field 
or by spin injection and thereby altering the output current, has so far been a peripheral 
topic. There are theoretical reports of spin filters and spin memories [§] in which the spin 
polarized current is controlled by the Zeeman splitting of the localized levels in the QD. An- 
other suggestion is a three-terminal system in which two of the leads are in anti-ferromagnetic 
order [Q. The source-drain current is manipulated by the magnetization direction in the 
third terminal. However, these studies are formulated in terms of single-electron properties 
and, in addition, they cannot be directly transformed into a time-dependent situation. To 
our knowledge, there is no theoretical report of inducing a large spin-polarization in the QD 
by simply spin-polarizing the conduction band. 

In this Letter we demonstrate, from a many-body approach, that there is a large spin- 
dependent renormalization of the levels in an interacting region, e.g. a QD, due to the 
magnetic properties of the leads, which could be used in magnetic sensor applications. By 
shifting the magnetic ordering in the leads, from paramagnetic to ferromagnetic, the levels 
in the interacting region experiences a spin split due to kinematic interactions with the 
conduction bands. In fact, as we will show, the conduction electrons with the spin projection 
a interact kinematically with the localized level of the opposite spin a. This effect, in turn. 
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causes a drastic increase (up to 150%) of the tunneling current through the interacting region 
when the conduction bands are changed from a ferromagnetic to paramagnetic ordering. 
Having this result at hand, we suggest a single electron device that is sensitive to the 
magnetic ordering in the contacts and is responding with an altered output current. 

To be specific, we are interested in an interacting region with a single level, that is 
taking part in the conduction, in the presence of a large Coulomb repulsion U. Such a 
system corresponds to the experimental reality of a QD in the Coulomb blockade state at 



low temperatures ksT < U, where /c^ is the Boltzmann contant, and small voltages pO 
The interacting region is coupled, via tunneling (mixing) interactions Vka, to two contact 
leads characterized by free electrons and the chemical potentials fi^ and fiR for the left (L) 
and the right (R) leads, respectively. A voltage applied over the system giving rise to a 
difference /xl — /ir 7^ results in a charge current from the higher to the lower chemical 



potential. The system can be realized with the degenerate Anderson Hamiltonian [|T^], with 
two conduction bands, in which the localized states in the interacting region are described 
hj Hd = J^pEpX^P, where the Hubbard operator X^^ = \p){p'\ represents the transition 
from the state \p') to \p) [0. The summation is taken over the state labels p e {0, t, |} 
(p = corresponds to the local vacuum whereas the doubly occupied state | t|) is excluded 
because of the large U). A conduction electron with the energy Eka- in the lead a = L, R is 
created (annihilated) by cj,^ (cka), ka G a. The Hamiltonian of the system can be written 
as 

n= e,A,c,, + Y.E,X^^ + Y.{v,A,X'^ + H.C.). (1) 

k(7eL,R P ka 

The dynamics of the operator X^'^ is given by the Heisenberg equation of motion, 

= [x°^ n] = A^oX^^^ + + x^cka + vi,x'-c,,] 

(note that a and a denote opposite spin projections). It is the last term in this expression 
that gives the mixing induced spin-dependent dressing of the localized level. For a clari- 
fication of this fact, let us consider the difference of the diagram expansions for standard 
Fermion operators and Hubbard operators. When one is dealing with the expectation value 
of operators, tin = (T^IliLiVi)^ Wick's decoupling of two operators r/i,?]^ results in the 
anticommutator {r/j,//^}. In the case of standard Fermion operators this anticommutator is 
a scalar, Cij, and the number of operators in the expectation value is, therefore, decreased 
by two, TT^v = J2cijFij7iN_2y where Fij = {TrjiTij) is a Fermion propagator. Now, in the case 
of Hubbard operators the anticommutator becomes yet again an operator, Pij. The number 
of operators in the expectation ttat value is, thereby, only reduced by one, ttn = Dij^^-i-i 
where Dij = {TrjiTij) is a propagator of Hubbard operators, implying that one has to make 
a decoupling also with P^j, since it remains in the expectation value vttv-i- The terms in 
the perturbation expansion coming from the decouplings with P^j give rise to the kinematic 
interactions and are characteristic for strongly correlated electron systems. In our particular 
case the decoupling in the first step gives {X^" ,J2s^'"cks} = (A:°° + X'''')cka + X^^Cka, 
whereas the kinematic interactions are, in the second step, generated by the commutator 
[X^", (X°° + X^^)cka + X^^Cko] = X^^Cka- Thls effect is clearly seen to arise solely due to 
correlations. 
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The density of electron states (DOS) for each spin projection a in the interacting re- 
gion is given by Po-(cc') = —\ j -aYmG^^^iuS), where the Green function (GF) Go-a'{t,t') = 
'-i){TSX^''{t)X'''%t'))/{TS), with S = exp{-i!^'^n'{t)&t). The Hamiltonian 



is a time-dependent disturbance to the system, by which a perturbation expansion of Go 
is generated through functional differentiation with respect to the fields U^{t) [0. The 
equation of motion for the GF of the interacting region is 



d_ 



+ [P..(t+) + i?..(t+)] / ' K(t,t")a.'(t",t')dt" 

JO 

+[p.^(t+) +i?.^(t+)] r\,{t,t")cit",t')dt". (2) 

Jo 

Here A°q = — Eq is the bare transition energy, AUaoif) = Ua(t) — Uoit) and 
R^^,{t) = 5^^,i5/5Uo{t) + i5/5U^:^{t). The expectation value P,,/ = (T{X°", X"'°})/(T5) = 
5„^iNq + N„i„ is the sum of the population numbers Nq and N„i„ corresponding to the 
transitions [00] and [cV], respectively. Physical quantities are obtained as U^{t) —>■ 0. 
In this limit, all expectation values which do not conserve the longitudinal component 
of the total spin vanish, although their functional derivatives may not. The propagator 
Va{t,t') = Y.keL,R I'^kcrl'^gkcrityt'), where gka is the GF of free electrons in the lead a. 

We look for a GF of the form [|l^ Gaa'it,t') = Daa'{t,t')Pcra'{i'), where the locator D^a' 
provides the essential physical information when all the P^ are approximated by constants. 
A more detailed study shows that taking into account effects of 6P^(t)/6U^'(t') only 



marginally modify our results. If we neglect all functional derivatives in Eq. (j^) the Hubbard- 
I approximation |T3| is recovered. By also calculating the first functional derivative of the 
GF, for which the only non- vanishing contribution is 

RMt^)G^^it",t') = -tD^^{f,t,)Toih,t2;t+)G^^ih,t'), 

we find the first order equation in the tunneling interaction V. Here, we have defined the 
zero vertex To{ti,t2;t+) = 6d^^{ti,t2)/SU^a{t+) = -6{ti -t2)6{t2-t+) |g, where d is the 
locator of the interacting region for vanishing tunneling interactions with the leads. The 
Dyson equation for the locator D„„ generated by the zero vertex, the loop correction, is 
graphically given by 



= — >- 




where the single and double straight lines symbolize the locators d and D, respectively. 
The wiggly line denotes the effective interaction V. We draw attention to the fact that the 
localized level with the spin projection a interacts kinematically only with the conduction 
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electrons of the opposite spin a and it is this effect that gives the possibihty of a large 
magneto resistance (MR). The renormalized transition energy Ag-o is given by the equation 

A AO _ V- L, |2 /(Aa0) - fjSkcr) 

where f{uj) is the Fermi-Dirac distribution function. Note that results similar to Eq. (Rf) have 



been obtained earlier by other methods for different models in equilibrium |15|. However, in 
none of these earlier studies the explicit spin-dependence on the conduction electrons of the 
opposite spin projection, which is present in Eq. @, was found. For constant mixing v and 
conduction band density of states the shift, given by Eq. @, clearly has a logarithmic 
divergence at the chemical potential of the lead a. Obviously then, for certain choices 
of parameters Eq. (^) has more than one solution, as illustrated in Fig. |l[ Since the 
renormalization of A^q (solid line) depends on the dressed transition energy Aj^q (dashed 
line), there may be several divergences around the chemical potential. All such solutions 
correspond to possible excitations of the QD. However, the interesting solution for each spin 
is that with the lowest energy. In Fig. ^ the dressed transition energies A^o (solid line) and 
Aj^o (dashed line) are plotted as a function of the spin polarization in the leads, defined as 
the fraction [W^ — W^)/W , where is the high energy cut off for the constant conduction 
band density of states in the lead a and W is half the bandwidth of the conduction band. 
Throughout this Letter we consider only the case when the polarizations in the two leads are 
the same. For non-polarized leads, the localized spin |- and spin |-levels collapse into a two- 
fold degenerate level. As the spin polarization in the leads becomes non-zero, the dressed 
transition energies for the two levels become distinct and as the polarization increases, the 
renormalization of the j-level decreases. In the limit of completely spin polarized conduction 
bands, the renormalization vanishes and A-fo —>■ A^g. 

The splitting of the localized levels in the interacting region due to the magnetic proper- 
ties in the conduction bands directly influences the tunneling current. Of particular interest 
is a comparison of the cases when the two leads are in either paramagnetic or ferromagnetic 
order, since these are the relevant states in magnetic sensors. Below we show that the dif- 
ferent magnetic phases of the leads imply a severe change in the magnitude of the tunneling 
current through the system. 

In the stationary regime, the tunneling current through the interacting region, symmet- 
rically coupled to the leads, is given by (for a detailed discussion see Refs. [p^,p3[1). 



J=tX\ r.(/L(^)-/i?(^))p.(^)de, 

where F,, = r^r^/(r^ + F^), F" = F^(^)|^„ = 2t^\v,\^pI and f^{e) = f{e - /ij, which has 
proven successful in the regime we consider |^^. In the given approximation the retarded 



GF is 

Pr 



era 



with A^o given by Eq. (g) and P^^ = Nq + N^, where = -l/n J f{u)lmGl^{iu)duj 
and Nq + N-^ + Nj^ = 1. The corresponding DOS is shown in Fig. |] where the spin | and 
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spin I are plotted on the positive and negative vertical axes, respectively. When the leads 
are in a paramagnetic state (dashed lines), the dressed transition energies coincide having 
an equal probability. For ferromagnetically ordered leads (solid lines) the A|o-transition 
becomes more likely (P^ = 0.69) than the Aj^Q-transition (P|| = 0.39). At the same time 
the transition to the |-level retains the strong influence from the spin | conduction electrons 
and therefore remains as large, or larger, as in the paramagnetic conflguration. 

As for the tunneling current through the system, there is a huge discrepancy in the 
current for a range of voltages in the two cases, displayed in Fig. ^ In Fig. ^ the current- 
voltage characteristics are shown for three cases in which the leads are in paramagnetic 
(dashed line) and ferromagnetic ordering with a minority spin percentage of 3.5 (solid line) 
and (dotted). For sufficiently small voltages the magnitude of the current is larger for 
ferromagnetic than for paramagnetic leads. As the voltage increases, though, the current 
becomes larger in the paramagnetic case and for certain voltages the change in the MR 
\Rfrn — Rpm\l Rpm cau be as large as 150%, a large number in view of existing experimental 
devices p5[| . 

The solid line in Fig. ^, describing a spin-dependent resonant tunneling behaviour, 
represents a situation where the spin polarization in the leads is such that there is only a 
tiny fraction of the minority spin state present. As the bias voltage is increased the bottom 
of that sub-band eventually separates from the corresponding level in the interacting region 
which gives a decreasing contribution from the minority spin to the tunnel current. This 
results in a tunnel current through the system that equals the current of the majority spin 
state only. The inset of Fig. ^ illustrates a non-intuitive and extreme case of this situation 
with very narrow conduction bands. Then, for a certain voltage range the current is decreased 
as the conduction bands are shifted from ferromagnetic to paramagnetic conflguration, giving 
an up to 45% inverse MR. 

In conclusion, using the Anderson model we predict that the localized level with a given 
spin state in a QD is strongly renormalized, via kinematic interactions, by the conduction 
band of the opposite spin state. For ferromagnetic leads the levels in the QD experience a 
spin split which results in a spin-dependent tunnel current through the system. We observe 
a change in the MR by up to 150% as the magnetic conflguration in the leads are changed 
from ferromagnetic to paramagnetic, suggesting that our flndings can be used in devicing 
magnetic sensors. The effect is non-trivial which is shown by the possibility of an, up to 
45%, inverse MR. 
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FIGURES 
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FIG. 1. The graphical solution of the renormalization Eq. (|3|) for the spin j (solid) and spin 
1 (dashed) level in the interacting region. The bare transition energy AJ^q = —0.1 relative to 
the chemical potential fia = 0, the coupling strength F" = 2Tr\va\'^p^ = 0.5 and the temperature 
kT = 0.175. The conduction band density of states = 1/2W = 1/100. The spin polarization in 
the conduction band is given by the lower cut offs — IF|* = — 50, — VFP = —12.5 (units: meV). 
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FIG. 2. The dressed transition energies A|o (solid) and A^o (dashed) as a function of the spin 
polarization (W?' — IFP) /W? in the conduction bands. The polarization in the two leads are equal. 
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FIG. 3. The equilibrium DOS of the interacting region for paramagnetically (dashed) and 
ferromagnetically (sohd) ordered leads. In the ferromagnetically ordered case there is an amount 
of 3.5% of the minority spin in the leads. 




FIG. 4. The tunneling current through the interacting region for paramagnetically (dashed) 
and ferromagnetically ordered leads. For the ferromagnetic ordering the minority spin percentage 
is 3.5 (sohd) and (dotted). The inset shows the situation where the conduction band width is 
2W = 2.8 meV. Iq = 4 + I^, = eF^/n. 
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